Spontaneous polygonization of multi-walled carbon nanotubes: perturbation analysis 
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Spontaneous polygonization for a multi- walled carbon nanotubes (MWCNTs) have been observed 
for about two decades. In present manuscript, this phenomenon is understood by the competition 
t-H ' between cohesion energy (with lattice mismatching effect) and curvature elastic energy of tubes 

t-H , inside the MWCNTs. Based on a continuum elastic model and perturbation analysis, the crucial 

. conditions for spontaneous polygonization of MWCNT is expressed by fundamental parameters of 

(N 1 MWCNTs, which are in reasonable agreement with all relevant experiments and computer simu- 

lations reported in literatures. Present approach can be used in further studies of crossectional 
r- ■ ■ geometry-dependent properties of MWCNTs. 

' PACS numbers: 61.48.-c, 61.46.Fg, 62.25.-g 
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\ — I. INTRODUCTION 

^ : 

Multi- walled carbon nanotubes (MWCNTs) arc concentric cylinders which can be thought of rolled from graphenc 
layers of graphites. The MWCNT was firstly discovered in 1991—^, after that, carbon nanotubes (CNTs) have aroused 
\ intense interest due to their extraordinary properties and promising potential for nanotechnology applications?. 

Because of CNTs' special atomic structures, they are expected to exhibit physical behaviors fundamentally different 
from those of graphite^. Electronically, a CNT can be either a semiconductor or a metal, depending on its diameter 
. and chirality, and the electronic band structure is sensitive to the cross-sectional geometry of the tube. Thus the 
shape transition of tube cross section may play an important role in the applications of nano-electronic devices and 
I , nano-electromechanical system (NEMS)^£. 

Some circular crosssectional tubes in a MWCNT can be polygonized in response to hydrostatic pressure^. The 
pressure induced polygonization is help to form CNTs bundles with 2D hexagonal closed-packed structures^. Similar 
pressure induced polygonization of CNTs has also been observed for CNTs under non-hydrostatic pressures^. Interest- 
ingly, polygonization can occur spontaneously for MWCNTs with large cross section even without radial pressure^—. 
Kiang et al£ used x-ray diffraction method, found polygonization (or graphitization) flattens the surface of the tube, 
decreases the interlayer space, as well as introduces sites with high curvatures. Wu et ali^- also studied the structural 
transition of MWCNTs by x-ray diffraction, suggested the larger the tube cross section is, the easier polygonization 
occurs. After heat treatment, they found critical diameter for polygonization being 50 nm (lower bound). Yoon et 

^ | alJ^ also studied the heat treated CNTs, showed that although normal CNTs are stable at low temperature, entropy 

due to Stone- Wales defects created during heat treatment make polygonal shape of larger CNTs more stable, and the 
radius of a polygonized tube is 50 — 100 nm. All the experiments, also at different conditions, suggested the shape 
transition of tubes is radius relevant, i.e., there is a critical radius for polygonization of tubes in a MWCNT. 

Spontaneous polygoniztion of tubes in MWCNTs has not been studied extensively by theorists yet, for that atomic 
simulations are time-consuming and continuum model usually lost some important information on lattice structure. 
The first continuum model for polygonization is proposed by Golovaty et al.—. They considered strong covalent bonds 
between the atoms in a graphene layer and weak bonds between the atoms in adjacent layers in their continuum 
theory of MWCNTs. Golovaty et al.'s work is based on the basic idea: the curvature-induced mismatching between 
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5h I the lattices of the adjacent graphene sheets can be reduced by flattening the walls at the expense of creating the 
. . i line defects that run parallel to the axis of the tube. Polygonization results from the competition between such two 
types of interactions, which is controlled by the diameter but not the number of walls in MWCNTs. Tibbetts et al*£ 
proposed a model considering many energy terms, such as bending strain energy, tensile strain energy, interfacial 
energy etc., showed that for a MWCNT with certain number of walls, optimal polygonal shape of CNT depends on 
the interlayer spacing and the chirality of the tubes. 

Although the physics of spontaneous polygonization is clear, all previous theoretical studies did not explicitly 
express the critical radius of this shape transition by basic physical parameters of CNTs. Golovaty et aFs model 
and geometric approach are beautiful but highly abstract, impossible to give quantitative results, which limits the 
application of their theory. Tibbetts et al. modeled the quasi-polygonal crosssectional CNTs as several flattened walls 
connected by small curved junctions at the corner of polygons, therefore their results rely on too many assumptions 
on the exact shape of the tube^. 

In present manuscript, we will give the explicit expression of the critical radius for polygonization using our geometric 
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theory of continuum clastic model. To use simple model illustrating basic physics, we start with a model containing 
elastic energy (curvature energy) and cohesion energy between adjacent layers, 

E — E cur + E co h- (1) 

For simplicity, wc consider a straight MWCNTs in which all tube cross sections are perpendicular to the axis. 
Polygonization is projected to a 2D planar geometric problem, as shown in Fig. [TJ A boundary of tube cross section 
is a closed plane curve described by the polar vector p = p(<fi), with 4> £ [0, 2ir) being polar angle. 

Following the approach in Ref. [l2|, we study the simplest case of two- walled CNTs, since polygonization is sensitive 
to the radius of tube, but not the number of walls. We assume that cross sections of these two tubes are concentric, 
similar to each other in the meaning of Euclidean plane geometry. Given the the closed curve for inner tube C\ being 
Pin = p{4>), obviously, p out — k p{4>) is the radius of the larger tube, and k is the geometric similarity ratio. The ratio 
for areas enclosed by C 2 and C\ is S 2 /Si = k 2 . 

Only curvature elastic energy is considered in present model, since the the strain energy corresponding to stretching 
of carbon-carbon bonds are much higher than curvature energy related to rotating the bonds, i.e., changing the bond 
angles. Curvature energy of a CNT tube have a form o f 14 i 15 , 

E c = y /(2-Ef) 2 dA + k j> KdA, (2) 

Here H and K are the mean curvature and Gaussian curvature for the tube surface. The bending constants k c = 
1.17 eV and k/k c = —1.56 was obtained by Ou-Yang et al*£, which can be expressed by basic bond parameters of 
Lenosky type lattice energy—, and in good agreement with the values reported by Tersoff et al*Z. The Eq. ([2]) can 
be projected to 2D elastic energy, 

E cuv = kJ i f k 2 ds (3) 

Here, k r is curvature of plane curve p((f>) as shown in Fig. [TJ s is arc parameter of the curve C\, L is tube length. 

In general, interlayer interaction can be obtained by the summation of inter-atomic potential between carbon atoms 
in adjacent layers^. For simplicity, we use the continuum form of interlayer interaction energy, which is proportional 
to the volume between inner and outer tubesi^, 

E coh = -goL (k 2 -l) <j> dA, (4) 

Here, dA is area element of area enclosed by curve Si. The —g = 7jAE c /d , AE C = —330 erg/cm 2 = — 2.04eV/nm 2 
is the interlayer cohesion energy of planar graphite obtained theoretically by Girifalco and La d 15 ' 18 , do ~ 0.34 nm is 
interlayer space of tubes in a MWCNT, which varies little with the increasing of number of shells ^ 11 ! 13 , thus we keep 
it being a constant. Two adjacent tubes may have different chiralitics, leading to incommensurate lattice structures. 
Typically, the interlayer interaction between two mismatched lattices is smaller than that for the commensurate ones. 
Therefore, the factor 77 is introduced to account for the mismatching effect between the interlayer lattices, < 77 < 1, 
and rj = 1 corresponds to the perfect matched lattices, i.e., graphite case. Exact value of mismatching factor rj should 
be determined numerically by first principle calculations. The total energy gives, 



E = E cuv + E coh = a j> k 2 r ds + p <j> dA, (5) 

Here, constants a, f3 and geometric similarity ratio can be expressed as 

k c L(k + l) . 2 . do 
— , f3 = -g L(k -1), k = l + — 



Equilibrium shape of the tube's cross section is derived from variation equation 6 J- = 0, and J- = E + A ds — Lq) . 
The perimeters of C\ are fixed due to the carbon-carbon bond lengths are unchanged during spontaneous polygoniza- 
tion. The restriction enters objective function of variation J- by the term with Lagrange factor A. Vanishing of SF 
gives the shape equation, 

akl(s) +2ak"(s) + /3 -2\k r (s) = 0, (6) 

where, "fc"(s)" denotes d 2 fc r /d.s 2 . 
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Obviously, circle is one of the solutions of Eq. ((6j) with k r = 1/Ri n , since this cubic algebraic cquaiton of 1/Ri n 
always has real roots. 

Polygonization means circular crosssectional tubes are instability, which we will discussed in detail by perturbation 
analysis. We consider a slight distortion of a circular cross-sectional tube pair with p((j>) = Rirr^, 



P (<j>)=Rin + Y, b ™ eim4 > 



(7) 



where <fi S [0, 27r), m = 0, ±1, ±2, . . . , ±oo, and b* n — b- m ensures p is real. High-order harmonic modes are presented 
in the expression of p((j>) as perturbation terms. The variations of the area enclosed by the curve Ci, length of 
circumference of C\ , and the curvature-elastic energy can be obtained as 1 ^, 



5 / d^ = irRi 
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The length of C\ is conserved, therefore, 5 J ds = 0, which provides one restriction on coefficient bo, 



op = _1 y 



Ri 



(11) 



Substitute the restriction Eq. (HTJ to Eq.([5|) and Eq. (fTTJj) . and the total energy shown as Eq. ([5]), gets, 



5 dA = nRf, 
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yields the deformation energy, 



5E ~ Triir 1 £ [2a (m 2 - l) - /?i??J ( m a _ j) 



m=0 



Rin 



The trivial case of m = 1, characterized by SE = 0, means a translation of the circle. The m-th harmonic 
deformation can happen when the coefficient of \b m /Ri n \ 2 in Eq. (|14[) becomes negative, which is, Ri n > R m , with 

1/3 

Rm = [2a(m 2 - l)/0\ ' , (m > 2). The occurring of m-order harmonic mode lowers the energy, as well as lowers 
the symmetry of tube from Coo to C m . Substitute the expressions of a and /3 to > R m gives, 



Rin ^ 



9o d 



x y m 2 — 1, 



(14) 



Experiments suggests hexagon is a typical quasi-polygonal shape of distorted CNTs, thus the critical radius for the 
occurring of hexagonal cross section is 



R 



critical 



i]AE c 



(15) 



Tubes with radius larger than R cr iticai ~ should be polygonized to hexagonal shape. For 77 = 0.1, R cr iticai = 

15 nm, which is in reasonable accord with the spontaneous polygonization phenomena of MWCNTs observed in recent 
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experiments^^. Based our calculation, wc confirme that mismatching of lattices (reflected by factor r\ < 1) plays an 
important role in polygonization, and mismatching effect can be amplified by the defects introduced in the process of 
heat treatment. 

In summary, we have studied the spontaneous polygonization of CNTs without applied hydrostatic pressure. The 
key step is the perturbation analysis for the instability of original circular crosssectional tube with large circumference. 
We only consider the bending energy coming from the interatomic interaction for atoms in the same layer and 
cohension energy due to interatomic interaction for atoms in adjacent layers, as well as the mismatching of lattices. 
With the length of circumference being fixed, the competition of bending energy and cohension energy lends to 
the polygonization. Wc obtained the the critical radius for tube occurring shape transition from circular shape to 
hexagonal shape, as shown in Eq. (|I5p . which is the first quantitatively expression for the critical radius of spontaneous 
polygonization, The results are in accord with recent experiments, and can be used to predict polygonization related 
phenomena. 
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FIG. 1: Cross sections of the simplest two-walled CNT. 



